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■ Abstract. Precise total cross-sections and invariant-mass distributions have been measured for photopro- 
duction of pion pairs off the proton producing pir'^Tv'^ and nvr'^'Tr" final states from the threshold region up 
to 800 MeV incident photon energy. Additionally, beam helicity asymmetries have been measured in the 
second resonance region (550 MeV - 820 MeV). The experiment was performed at the tagged photon beam 

k>( , of the Mainz MAMI accelerator with the Crystal Ball and TAPS detectors combined to give an almost 

■ i-jT solid-angle electromagnetic calorimeter. The results are much more precise than any previous measure- 
' ments and confirm the chiral perturbation theory predictions for the threshold behavior of these reactions. 

In the second resonance region, the invariant-mass distributions of meson-meson and meson-nucleon pairs 
are in reasonable agreement with model predictions, but none of the models reproduce the asymmetries 
for the mixed-charge channel. 

PACS. 13.60.Le Meson production - 14.20.Gk Baryon resonances with S=0 - 25.20.Lj Photoproduction 
reactions 



1 Introduction 

The photoproduction of pion pairs has attracted a lot of 
interest over recent years in view of different properties of 
the strong interaction. In particular, it can contribute to 
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sensitive tests of chiral perturbation theory, to the inves- 
tigation of non-ground-state decays of nucleon resonances 
via reaction chains like R — > R'tt Nirn {R, R' nucleon 
resonances, N = n^p), and to the analysis of in-medium 
properties of hadrons like the tr-meson. The first topic is 
related to the threshold behavior of double pion produc- 
tion off the nucleon, which has been predicted by chiral 
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perturbation theory P^I14) . The second is important for 
the extraction of nucleon resonance properties from the 
electromagnetic excitations of the proton and the neu- 
tron. Finally, the interpretation of modifications of the 
pion-pion invariant-mass distributions in nuclear matter 
requires input from the elementary reactions off the free 
nucleon. 

Most of the recent progress in the theoretical treatment 
of the strong interaction comes from two different lines of 
research: (1) the development of the numerical methods 
of lattice gauge theory, approaching the non-perturbative 
regime downwards from large quark masses (see e.g. [1] 
[21[3]) and (2) the methods of chiral perturbation theory 
(ChPT) [31[S] , extrapolating upwards from small momenta 
and vanishing quark masses. The latter, originally based 
on the approximate Goldstone boson nature of the pion 
and later extended to the nucleon sector [Bl|7] , had a large 
success in correctly describing the threshold behavior of 
photon-induced 7r° production off the proton. Agreement 
with experiment [8ll9lfT0] was achieved due to the con- 
tribution of pion loop diagrams, which were not consid- 
ered in the derivation of the older low-energy theorems 
(LET). Subsequently, different observables in single pion 
photoproduction have been measured and interpreted in 
the framework of chiral perturbation theory. Examples 
are the unitary cusp in the s-wave i?o+ multipole of tt*^- 
photoproduction at the tt"*" production threshold fll' and 
the polarizability of the charged pion [12, studied via the 
IP ~^ "fTT~^n reaction. In the meantime, it was discov- 
ered [T^fTl] that the chiral loop diagrams play an even 
more important role for double pion production, opening 
a new door for ChPT tests. An unexpected prediction was 
that, very close to threshold, the double 7r° channel is so 
much enhanced by loop corrections that its cross section 
rises above the charged and mixed-charged channels. This 
is completely counterintuitive. At higher incident photon 
energies the reactions involving pairs with charged pions 
(Tr^Tr"*", 7r~7r+) have much larger cross sections, and in sin- 
gle pion production in the threshold region, the pTr" final 
state is strongly suppressed with respect to mr^ . How- 
ever, in double pion production, although contributions 
of order one vanish for neutral pion pairs, the corrections 
of order Afjr are predicted to be much larger than for the 
pairs with charged pions \2A, ■ This result is stable when all 
next-to-leading order terms of order are considered. 
At this order also resonance contributions come into play 
via the N*Nnn s-wave vertex. The largest contribution 
by far involves the Pii(1440) resonance. Taking all contri- 
butions together, the leading-order chiral loop diagrams 
contribute roughly 2/3 of the total 27r° yield, making this 
channel the ideal testing ground for these contributions, 
which in most other reactions account only for small cor- 
rections. The calculations [14] predict the following thresh- 
old behavior for 7p — >■ ptt^tt": 



^tot(^7) =0-6 nb- 
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where Ej is the incident photon energy in MeV and E^^^^ = 
308.8 MeV is the production threshold. The largest uncer- 
tainty in this prediction comes from the s-wave coupling of 
the Pii(1440) resonance, which was taken from an analy- 
sis of the ttN ttttN reaction [TS] . Equation ([1]) uses the 
central value of this coupling. If instead the upper limit is 
used, the coefficient in the equation rises from 0.6 nb to 
0.9 nb. 

The threshold behavior of double tt*^ production was 
also investigated in the framework of the Gomez Tejedor- 
Oset (Valencia) model . Here, it was found p/7j that the 
low-energy cross section is considerably enhanced when 
final-state interaction (FSI) contributions are taken into 
account in addition to the pure tree-level treatment in 
|16) . In particular the re-scattering of 7r+7r^-pairs pro- 
duced via the Zi-KroU-Ruderman diagram into tt'^tt'^ gives 
a considerable contribution. Altogether, the FSI effects al- 
most double the threshold cross section for the 27r'' chan- 
nel, which can be regarded as a kind of remnant of the 
chiral loop effects discussed above. However, very close to 
threshold, their final result is still smaller than the values 
predicted by the chiral theory. 

At higher incident photon energies, double pion pro- 
duction is of interest in view of the decay of nucleon res- 
onances. It is a well-known problem that far fewer nu- 
cleon resonances have been observed in experiments than 
are predicted by quark models. A possible explanation 
could be the experimental bias against states that cou- 
ple only weakly to Ntt, introduced into the data base 
by the dominance of elastic pion scattering experiments. 
Photon-induced reactions can remove this bias for the ini- 
tial state, but then final states other than Ntt must be 
investigated. As discussed in detail in [TH] three alterna- 
tive decay modes of nucleon resonances can be studied via 
double pion production: the sequential decay via an inter- 
mediate excited nucleon state (for all types of pion pairs) , 
the decay into the Np channel (for tt'-'tt^ or tt+tt") and 
the emission of the cr- meson (for tt'^tt'-' or tt+tt^). The reac- 
tion 7p — p7r°7r" has been recently investigated in detail 
in view of such resonance contributions with the Bonn- 
Gatchina (BoGa) coupled-channel model up to incident 
photon energies of 1.3 GeV [l9l[20] . 

Much more data are available at lower incident photon 
energies in the second resonance region comprised of the 
Pii(1440), 511(1535), and L>i3(1520) states. Total cross 
sections and invariant-mass distributions of the tttt- and 
the TrA^-pairs have been measured with the DAPHNE and 
TAPS detectors at the MAMI accelerator in Mainz pTl 
[22l[23l[24l[25l[26l[27l[28] . at GRAAL in Grenoble (also hn- 
early polarized beam asymmetry) [2^130) . and with the 
CLAS detector at JLab (electroproduction) fHT] . Recently 
polarization observables were also measured at the MAMI 
accelerator [321|31IM1I351I3S] and at the CLAS facility [37]. 
However, it is somewhat surprising that, despite all these 
efforts, the interpretation of the data is still controversial 
even at low excitation energies, where only a few well- 
known resonances contribute. All models agree that the 
production of charged pions involves larger background 
contributions from non-resonant terms like pion-pole dia- 
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grams or terms of the zi-KroU Rudermann type. The ex- 
treme cases are the tt^tt^ final state, which at moderate 
incident photon energies is dominated by such contribu- 
tions, and the tt^tt^ final state with only small background 
contributions. The latter is thus better suited for the study 
of resonance contributions. 

However, the reaction models do not even agree for the 
dominant contributions to double tt'^ production. In the 
model of the Valencia group [16,38,39 the most important 
contribution comes from the Z3i3(1520) — > Att^ pn'^w'^ 
reaction chain. In the model of Laget and co-workers 
a much more prominent role is played by the Pii(1440) — )• 
Na decay. The recent Bonn-Gatchina analysis finds a large 
contribution from the £'33(1700) state, which is almost 
negligible in the other models. The effective Lagrangian 
model from Fix and Arenhovel (Two-Pion-MAID [3D]) is 
also dominated by the D13 resonance, but strongly under- 
estimates the total cross section at incident photon ener- 
gies below 700 MeV. The results for the helicity depen- 
dence of the cross sections obtained from the experiments 
on the Gerasimov-Drell-Hearn sum rule [5^I33[IM] show a 
dominance of the (J3/2 component. This is in line with the 
excitation of D13 or D33 resonances and limits possible 
Pii contributions. 

The situation is even more complicated for the mix- 
charged ri7r°7r+ and p7r°7r~ final states. Early model cal- 
culations in the framework of the Valencia or the Laget 
model (see e.g. [16]) could not even reproduce the total 
cross section. First experimental data for the invariant- 
mass distributions of the pion pairs [^[?7| indicated an 
enhancement at large invariant masses, which was inter- 
preted as a possible contribution from p-mesons, e.g. via 
the I?i3(1520) — > Np decay. Introduction of this pro- 
cess into the models improved substantially the agreement 
with experiment [SS l lM l liD ] . 

The behavior and interpretation of the low-energy dou- 
ble pion production reactions are also relevant for differ- 
ent aspects of the much discussed in-medium properties of 
hadrons. The above mentioned coupling of the Di3(1520) 
resonance to the Np channel has been discussed in [57] 
as a possible explanation of the strong suppression of the 
second resonance bump in total photoabsorption by nuclei 
[41] . For the free nucleon this structure is dominated by 
this Di3 resonance. If it couples strongly to Np and the 
p is broadened or shifted to lower mass in nuclear matter, 
then also the nucleon resonance will be broadened and the 
bump structure will be suppressed for nuclei. The study 
of the line-shape of this resonance via single 7r° photo- 
production off nuclei could not establish such a nuclear 
broadening [32] ■ However, the interpretation of the data 
is greatly complicated by final-state-interaction effects. 

The other important in-medium effect discussed in the 
context of low-energy double pion production is the be- 
havior of the (T-meson in nuclear matter. The masses of 
the ~ 0^ pion and its chiral partner the = 0+ 
(T-meson are very different in vacuum, which is a well- 
known manifestation of chiral symmetry breaking. The 
masses should become degenerate in the chiral limit, so 
that one naively expects a density dependence of the mass 



due to partial chiral restoration effects. Results 17 ob- 
tained in the framework of the Nambu- Jona-Lasinio model 
indicated indeed a strong drop of the a mass as a func- 
tion of nuclear density which is already significant at nor- 
mal nuclear density po at which the pion mass remains 
stable. The predicted effect for the production of meson 
pairs off nuclei is a downward shift of the invariant-mass 
distributions of scalar isoscalar meson pairs. This predic- 
tion has been experimentally investigated with pion- and 
photon-induced double pion production reactions 1441 1451 
136,47,48,49,50,51 . The CHAOS collaboration i44,'451|3g 
H71H5] reported such a downward shift for tt+tt" pairs 
with respect to tt^tt"*" pairs from pion induced reactions. 
The Crystal Ball collaboration at BNL [3^ observed a 
low-mass enhancement of strength for heavy nuclei in 7r~- 
induced 7r°7r° production. In photon-induced reactions, a 
downward shift of the invariant-mass distributions of tt^tt^ 
pairs with respect to tt^tt^ pairs has been measured by the 
TAPS collaboration [501151] . However, intricate final-state- 
interaction effects [5T] complicate the interpretation of the 
results. They require detailed studies in the framework of 
models, which have to rely on precise input for the ele- 
mentary cross sections off the free nucleon, in particular 
at low incident photon energies. 

The present paper reports a simultaneous, precise mea- 
surement of the total cross section and the invariant mass 
distributions of the 7p — )■ ptt^tt" and 7p — >■ titt'^tt^ reac- 
tions from as close as possible to threshold up to the sec- 
ond resonance region. In addition, the use of a circularly 
polarized photon beam allowed the measurement of the 
beam-helicity asymmetry as a function of the azimuthal 
angle between the reaction and production planes. The 
results for this asymmetry, which are very sensitive to de- 
tails of the reaction models, have already been published 
in a preceding Letter [36| . 

The paper is organized as follows. A brief description 
of the experimental setup is given in Sec. [2] The data 
analysis, including calibration procedures, particle identi- 
fication, Monte Carlo simulations of the detector response, 
and estimates of systematic uncertainties, is summarized 
in Sec. [3] In Sec. |3] (Results) first the threshold region is 
discussed, particularly in view of chiral perturbation the- 
ory predictions, followed by a comparison of the results 
from several reaction models to the observables at higher 
incident photon energies. 

2 Experimental setup 

The experiment was performed at the Mainzer Mikrotron 
(MAMI B) accelerator [5^155] using the Glasgow tagging 
spectrometer 54','55' and an almost Att covering electro- 
magnetic calorimeter combining the TAPS [571155] and 
Crystal Ball (CB) [5S] detectors. Here we give only a short 
summary of the main parameters of the setup. Details can 
be found in [HI], which used the same data set for the in- 
vestigation of the 7p — > pir'^y reaction. 

A schematic drawing of the main components of the 
detector system is shown in Fig. [T] The data were taken 
with a 4.8 cm long liquid hydrogen target, which was 
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mounted from the upstream side in the center of the Crys- 
tal BalL The Bah is composed of 672 Nal crystals covering 
the full azimuthal range for polar angles between 20° and 
160°. A Particle Identification Detector (PID) [35] (24 
plastic scintillator bars of 31 cm length, 13 mm width and 
2 mm thickness) and two cylindrical multiple wire propor- 
tional chambers (MWPCs) [60 were mounted inside the 
Ball in cylindrical geometry around the target, covering 
the same polar angle range. The forward angular range 
from 20° down to 1° was covered by the TAPS detector 
[571155] with 510 BaF2 crystals arranged as a hexagonal 
wall and placed 1.75 m downstream from the target. Each 
module of this detector had an individual plastic scintilla- 
tor (5 mm thickness) of the same hexagonal geometry in 
front for the identification of charged particles. 




Fig. 1. Experimental setup of the electromagnetic calorimeter 
consisting of TAPS and Crystal Ball detector. Detectors for 
charged particle identification are mounted inside the Crystal 
Ball around the target (PID and MWPCs) and in front of the 
TAPS forward wall (TAPS Veto-detector). 

The energy resolution for electromagnetic showers of 
both CB and TAPS is approximately given by [5S1I55] 

OE ^ 2-3% 
E ^ ^E/GeV 

A more precise parameterization for TAPS can be found 
in [551. The angular resolution was better than 1° FWHM 
for TAPS and FWHMe = 4.5° - 7° for polar and FWHM^ 
= FWHMe/sin(0) for azimuthal angles in the CB. 

The photon beam with energies up to 820 MeV was 
produced by bremsstrahlung from the 883 MeV electron 
beam of the MAMI accelerator. The energy of the incident 
photons was determined event-by-event by the Glasgow 
photon tagger [541155] with an energy resolution of approx- 
imately 2 MeV full width. The electron beam was longi- 
tudinally polarized with a polarization degree of (82±5)% 
so that the photon beam carried a circular polarization 
determined by the photon-energy dependent polarization 
transfer factor [62] . The incident photon flux was derived 
from the number of deflected electrons per tagger chan- 
nel, which were counted with live-time gated scalers. The 
tagging efhciency, i.e. the fraction of correlated photons 



that pass through the collimator and impinge on the tar- 
get, was periodically determined with special tagging ef- 
ficiency runs. This was done by measuring directly the 
photon beam intensity with a total absorbing lead-glass 
counter, which was moved into the photon beam at re- 
duced intensity. During the normal data-taking runs the 
photon beam intensity was monitored in arbitrary units 
with an ionization chamber at the end of the beam line. 
After normalization to the tagging efficiency runs these 
data were used to correct the time dependence of the tag- 
ging efficiency, which however was quite stable (varying 
between 30% - 34%). 



3 Data analysis 

3.1 Particle identification and reconstruction 

In general, electromagnetic showers produce signals in ex- 
tended 'clusters' of detector modules in the calorimeters. 
The first step of the analysis therefore combined all hits of 
adjacent crystals into 'clusters' and determined their en- 
ergy sums and their energy weighted geometrical centers 
of gravity. In the next step these 'clusters' were assigned to 
different particle types with the methods discussed below. 



3.1.1 TAPS forward wall with veto-detector 

The response of the TAPS-detector to electromagnetic 
showers originating from photons is discussed in detail 
in (55). The separation between photons, charged pious, 
and recoil nucleons in TAPS is based on three methods. 
The plastic scintillators from the veto-detector distinguish 
between charged and neutral hits. A hit was assigned as 
charged when the veto of any module in the cluster or the 
veto of any neighbor of the central module (module with 
highest energy deposit) had fired (even if the neighbor 
module itself had no signal above threshold). The latter 
condition applies to cases where a charged particle with 
large impact angle passed through the edge of a veto and 
then deposited its energy in a neighboring BaF2 module. 

Photons and recoil nucleons can be separated by a 
pulse-shape analysis (PSA) based on the scintillation prop- 
erties of BaF2. The crystals emit scintillation light at two 
different wavelengths with very different decay constants 
(t/ = 0.76 ns, Ts — 620 ns) and the intensity ratio of 
the two components depends on the nature of the inci- 
dent radiation. This is routinely explored by integrating 
the output signals over a short (w50 ns) and a long gate 
(w2 fis) period. The two energy signals were calibrated in 
a way that in a plot of long versus short-gate signal pho- 
tons appear along the 45° line. Since the fast component 
is suppressed for recoil nucleons, they appear at smaller 
angles (see insert in Fig. [5|). For practical purposes the 
signals were parameterized in polar coordinates R, (f> via 

R=^{E^+El), = tan-i(^|i) , (3) 
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Fig. 2. Pulse-shape analysis for the TAPS detector. The insert 
shows the signal integrated over a short time gate versus the 
total signal. The main plot shows the pulse-shape signal in 
polar coordinates. Solid (blue) line: position of photon peak, 
dashed (red) line: 3a limit cut between photons and nucleons. 

where Es, Ei are the short and long gate energy signals. 
The cut for the identification of photons was then defined 
in small slices of the spectra in the radius R projected 
onto the (/)-axis. These spectra were fitted by Gaussian 
peaks plus a first order background polynomial (see ^63^ 
for details). The solid (blue line) in Fig. [2] indicates the 
position of the Gaussian peaks and the dashed (red) line 
the 3(7 limit. Entries at the right hand side of the dashed 
line were accepted as photons. 
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50 100 150 200 250 300 350 400 450 500 

Energy [MeV] 

Fig. 3. Time-of-flight versus energy analysis for nucleons de- 
tected in TAPS. The offset of the time calibration is such that 
photons appear at zero time-of-flight. 

A further possibility to separate protons and photons 
in TAPS is a time-of-fiight versus energy analysis. This is 
demonstrated in Fig. [31 The selected events had exactly 
two photons in the Crystal Ball with an invariant mass 
close to the 7r° mass and one further hit in TAPS. This 
means that, apart from small background contributions, 
the reaction 7p — pii^ with the recoil proton in TAPS was 



selected. In the figure the difference between the average 
time of the two photons and the time-of-fiight for the hit 
in TAPS is plotted versus the energy deposited by the hit 
in TAPS. Proton hits are confined in a well-defined curved 
zone. The indicated limits define which hits are assigned 
as protons. The PSA and time-of-flight versus energy anal- 
ysis are complementary since the first is most efficient for 
large proton energies, while the latter is optimal for the 
smallest energies. In principle, it is possible to identify 
with these tools photons, neutrons, protons, and charged 
pions (see [ST | for details, the pions form a separate band 
in the time-of-flight versus energy spectra). However, as 
discussed in Sec. 13.21 for the present analysis they were 
only used to identify a very clean sample of photons in 
TAPS. 



3.1.2 Crystal Ball with PID and MWPCs 

Charged particles hitting the CB must traverse the PID 
and the MWPCs. The PID was used to identify protons 
and charged pions. This was done by a AE — E analysis 
that compared the differential energy loss of the particles 
in the PID to the total deposited energy in the CB. A 
typical spectrum is shown in Fig.S] The bands for protons 
and charged pions are clearly separated. 



1 1 1 1 1 1 1 1 1 1 1 1 1 




- Electrons 


Protons i 






-. T. . i' 1 1 1 1 1 1 'i 1 



"0 50 100 150 200 250 300 350 400 450 500 

Crystal Ball Energy [MeV] 

Fig. 4. AE — iJ-analysis using CB and PID. The energy de- 
position in the PID is plotted versus that in the CB. 

The energy deposition of hadrons spreads over fewer 
detector modules than electromagnetic showers and in 
many cases the energy deposit is confined to a single crys- 
tal. This means that the angular resolution for hadrons 
(w 10° for O) is worse than for photons for which the 
center of gravity of the extended cluster defines the im- 
pact point better than the detector granularity. Therefore, 
for charged particles the angular information from the CB 
was replaced by the tracking information delivered by the 
MWPCs whenever such information was available (if not 
the reconstruction from the CB clusters was used). Using 
the intersection points of the particle trajectory with the 
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two MWPCs an angular resolution of FWHMe = 3° - 
5.5° and FWHM^ = 3.3° was obtained. The angular reso- 
lution and the detection efficiency of the MWPCs was ex- 
perimentally determined with penetrating cosmic muons 
and the reactions 7p — >■ p7r° and 7p — t- n7r+. For the lat- 
ter two, events with protons, respectively charged pions, 
identified by the E — AE analysis were selected and it was 
measured for which fraction of those events the MWPCs 
had responded. The total efficiency of the MWPCs deter- 
mined this way was 90% for protons 79% for 7r+ mesons. 
Note, however, that this efficiency does not directly enter 
into the total detection efficiency since tt^ mesons with- 
out MWPC information were not discarded but had only 
a moderately lower angular resolution form the CB cluster 
reconstruction (see section for details of the efficiency 
simulations) . 



3.2 Reaction identification 



3.2.1 The reaction 



For this reaction events with exactly four neutral hits and 
one or no candidate for the recoil proton were accepted. 
Photons and proton (if detected) had to fulfill the above 
discussed identification criteria. Detection of the recoil 
protons was not required for two reasons. Protons from 
reactions at low incident photon energies have low kinetic 
energies and were mostly stopped in the target or other 
material. Requiring proton hits would have drastically re- 
duced the overall detection efficiency and would have not 
allowed a measurement of the reaction close to threshold. 
At higher incident photon energies, protons were detected 
with good detection efficiencies. However, when detection 
of the proton is required, the simulation of the detection 
efficiency becomes more involved and more model depen- 
dent. This would have unnecessarily introduced an addi- 
tional systematic uncertainty. 
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Fig. 5. Invariant masses of the 'best' combination of four pho- 
tons to two pairs for the 7p 



pTT^Tr" reaction. 



The accepted events were subjected to a combined in- 
variant and missing mass analysis. In the first step, the 
four photons were combined into the three possible dis- 
junct pairs. For each of the combinations the invariant 
mass of the two pairs was calculated and the 'best' com- 
bination was chosen with a test minimizing 



2 

k=l 



(m^-y(fc) — niT^aY 



(4) 



where m^o is the pion mass and the the invariant 

masses of the photon pairs with their uncertainties Anijy . 
The two-dimensional spectrum of the invariant masses of 
the best combinations in Fig. [5] shows a clear peak at the 
position of the tt" invariant mass. The small background 
under the peak was determined from a side-bin analysis 
and subtracted. 

Subsequently, the nominal mass of the pion was used to 
improve the resolution. Since the angular resolution of the 
detector system is much better than the energy resolution 
this was simply done by replacing the measured photon 
energies by 



E'12 = Ei^ 



(5) 



where Ei^2 are the measured photon energies, E'^ ^ the re- 
calculated energies, •nij.o is the nominal tk^ mass, and m^.^ 
the measured invariant mass. 
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Fig. 6. Missing mass spectra for the 7J3 — >■ pjr tt reaction in 
the threshold range and at high incident photon energies. The 
solid (red) histograms are GEANT simulations of the detector 
response. Vertical lines indicate the applied cuts. 



In the final step of the analysis the mass of the recoil 
proton, which was treated as missing particle regardless 
of whether it was detected or not, was compared to the 
missing mass of the reaction using 



Z\m(7r7r) 



P. 



Pp - P^o - P^ 



(6) 



where rup is the nucleon mass, P^, Pp, P^o ^ are the four- 
momenta of the incident photon, the initial-state proton 
(which was at rest), and the produced 7r°-mesons. The 
results for the two most critical energy regions - close 
to threshold and at highest incident-photon energy - are 
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shown in Fig. El The critical point in the threshold region 
is the very low cross sections of double 7r°-production. 
However, the spectrum shows a rather clean peak. At 
higher incident photon energies, earlier experiments [221 
US] were plagued by background from the r] — )• Stt" decay. 
The large solid angle coverage of the present experiment 
results in a negligible probability to lose two out of six 
photons so that this background did not contribute sig- 
nificantly. Only a tiny remnant signal is visible around 
missing masses of +100 MeV. 

In the final step, the kinetic energy and momentum of 
the recoil proton (no matter if detected or not) were recon- 
structed via the overdetermined reaction kinematics from 
the incident photon energy and the measured momenta of 
the two pions. 

3.2.2 The reaction — > rni^TT^ 

For this reaction events with two or three neutral hits and 
a charged hit in the CB were selected. Detection of the re- 
coil neutron was allowed but not required. The charged hit 
in the CB had to pass the AE — E analysis as a charged 
pion. Charged pions in TAPS were not considered. The 
reason is that the identification via the AE — E analy- 
sis using the PID was cleaner than the identification in 
TAPS where only time-of-fiight versus energy could be 
used to distinguish charged pions from protons (charged 
pions would appear in Fig. [31 between protons and pho- 
tons; see [5T] for details). Since protons misidentified as 
charged pions were the most important background source 
(see below) , events with cliarged-pion candidates in TAPS 
were discarded. This introduced of course an additional 
systematic uncertainty into the simulation of the detec- 
tion efficiency since a small part of the solid angle was not 
covered. 

In the first step of the analysis the invariant mass of 
the photon pair was calculated. When three neutral hits 
had been detected in the CB, where neutrons and photons 
cannot be distinguished, again the best combination was 
chosen. The resulting invariant mass spectrum was very 
clean and a cut was applied for invariant masses between 
115 MeV and 160 MeV. As for the double 7r° channel the 
nominal mass of the tt*' was then used to re-calculate the 
photon energies from Eq. ([5]). 

Treating again the recoil nucleon as a missing parti- 
cle, the missing mass can be calculated from Eq. (jG]), re- 
placing one of the 7r*'-mesons by the 7r+ and the final- 
state nucleon by the neutron. The result is shown in the 
center row of Fig. [7l The peaks around zero correspond 
to the 7p n7r°7r+ reaction. The background level is 
quite high, especially in the threshold region. The prob- 
lem arises from protons from the 7p pi:^ reaction that 
leak in the AE — E analysis (see Fig. U]) into the 7r+ re- 
gion. The probability for this leakage is small, however 
at energies below 400 MeV the cross section for single tt" 
photoproduction is larger than for the 7p nTr^Tr"*" re- 
action by roughly three orders of magnitude. In order to 
reduce this background, first the missing mass was calcu- 
lated under the hypothesis of single tt'^ photoproduction; 



i.e. the tt^ candidate was assumed to be a recoil proton 
and the missing mass 

Am{TT) = \P^ + Pp-P^o\~mp , (7) 

was formed. The result is shown in the upper row of Fig. [71 
In this spectrum, the background from single tt*^ produc- 
tion sits in the peaks around zero, while events from 7p ^• 
TiTT^TT"*" appear at large missing masses. For further analy- 
sis only the events in the shaded areas were accepted. The 
bottom row of Fig. [3 shows the two-pion missing mass 
after this cut on the one-pion missing mass had been ap- 
plied. These signals were practically background free with 
only a tiny contribution from 7] tt'^tt+tt" appearing at 
high incident photon energies and large missing mass. As 
for the double tt*^ case the measured distributions agreed 
very well with the simulated line shapes. 




Am(7i7i)[MeV] 



Fig. 7. Upper row: missing mass calculated from the tt kine- 
matics for the hypothesis of the 7p — > ptt" reaction. The peaks 
centered around zero are related to background from this reac- 
tion. The shaded areas were selected for further analysis of the 
7P nTT^TT^ reaction. Center and bottom row: missing mass 
calculated from vr" and tt"*" kinematics for the hypothesis of the 
7P — >■ nTT^TT^ reaction. Center: all events, bottom: after cut on 
shaded areas in Am{'K). Solid (red) histograms: simulation of 
detector response. Vertical lines: applied cuts. 

At low incident photon energies there was another small, 
but significant background component that could not be 
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removed completely with the invariant-mass analysis. It 
arises from the 7p — ^ ptt^tt" reaction. Very slow tt" mesons 
can be stopped inside the liquid hydrogen target and then 
be captured by a proton to form pionic hydrogen. The pi- 
onic hydrogen can subsequently decay by charge exchange 
via the Tr~p — >■ mr'^ reaction. Since this involves 7r~ that 
started with very low momenta, their energy loss is not 
significant within the experimental resolution and they 
cannot be discriminated by the missing mass analysis. 
However, since the secondary tt*^ mesons decay practically 
(within experimental resolution) at rest in the laboratory 
they can be eliminated in the pion kinetic energy spectra 
(peak at zero energy) or even in the opening angle spectra 
of the pion decay photons (peak at 180°). 

In the final step again the kinetic energy and momen- 
tum of the recoil neutron were reconstructed from the re- 
action kinematics, regardless whether the neutron was de- 
tected or not. 



3.3 Extraction of cross sections and systematic 
uncertainties 

The absolute normalization of the cross sections was based 
on the measurement of the incident photon flux, the target 
density, the two-photon-decay branching ratio of the tt*^ 
which is (98.823±0.034)% [M], and the efficiency of the 
detector system. 

The photon flux was determined as explained in Sec- 
tion [21 However, for the present experiment an unsolved 
problem occurred in the read-out electronics, which led to 
a staggered pattern in all types of Nd{E^)/N^{E^) ratios, 
where Nd are the counts for any reaction observed in the 
calorimeter and is the photon flux. The pattern oc- 
curred in groups of four tagger channels with a maximum 
amplitude of ±3.6%. It was reduced to the average val- 
ues by applying correction factors and a total systematic 
uncertainty of 5% was adopted. 

The target surface density was 0.201 nuclei/barn with 
a systematic uncertainty of 2%. Contributions from the 
target windows (in total 125 /im Kapton) were determined 
with empty target measurements and subtracted. For the 
first part of the beam time, a build-up of ice on the down- 
stream target window due to water permeation through 
the outer target tube was observed. The relative thickness 
of this layer was monitored using reactions with protons 
at large polar angles (> 80°), which can only arise from 
the heavy nuclei in the target windows and the ice but 
not from reactions on the liquid hydrogen. The track re- 
construction with the MWPCs was used to measure the 
intensity of such events from the different windows. Sub- 
sequently, the results from a measurement with a water 
target, normalized to the thickness of the ice layers, were 
used to subtract this background. For the second half of 
the beam time, this problem was eliminated by a modi- 
fication of the target. The correction for the 7r°7r° cross 
section due to the ice layer for the first part of the beam 
time amounted to « 10%. We estimate the total system- 
atic uncertainty due to this correction is below the 2% 
level. 



The systematic uncertainty due to the elimination of 
background via the invariant- and missing-mass analy- 
sis (including the agreement between the observed line 
shapes and simulations) is estimated to be about the ±3% 
level for the tt^tt" final state and at ±7% for the 7r°7r+ 
final state. It is larger for the latter due to the addi- 
tional background from stopped pious undergoing charge 
exchange and because there is no invariant mass filter for 
the charged pion. 

The detection efficiency of the experimental setup was 
simulated with the GEANT3 program package [65]. All 
details of the experimental setup (active detector compo- 
nents as well as support structures and other passive ma- 
terials) were included in the simulations and the results 
were tested in detail against known cross section data. The 
detector response to electromagnetic showers induced by 
photons is very well known. The systematic uncertainty 
for single photon detection is so small, that even at very 
high statistical accuracy no significant deviations of cross 
section data constructed from the 77 Stt" 67 decay 
from the world data base was observed [66 . This was ex- 
ploited to increase the statistical quality of //-production 
data by simultaneous measurements of the 27- 67-decay 
modes (see e.g. [SZ])) again without observation of any 
systematic differences. For the charged pious, in addition 
the MWPCs must be considered, which are not routinely 
included in the GEANT simulation. They were treated in 
the following way. Efficiency and angular resolution of the 
chambers were experimentally determined (see sec. 13.1.2)) . 
In the simulation, position information from the chambers 
was generated according to these experimental parameters 
from the, in the simulation exactly known, tracks of the 
charged pious. Whenever such information was available 
it replaced in the analysis, like for the measured data, the 
position information from the CB cluster. If not, the MW- 
PCs were ignored. For the detection efficiency of charged 
pious, this is only a second order effect. The efficiency is 
mainly determined by the GEANT3 simulation of PID 
and CB. The quality of GEANT3 simulations for charged 
pious in a calorimeter has been for example studied in |51) , 
using the rj — > tt'^tt+tt" decay. Also the agreement between 
the line-shapes of measured and simulated data (cf. Fig. 
[7]) demonstrates the high quality of the simulations. 

The dominant uncertainty of the detection efficiency 
is related to the event generators used for the simulation, 
which should reflect the 'true' kinematic correlations be- 
tween the two pious. This is discussed in detail below. 
For the 7p — > pir'^'K^ reaction, two different event gen- 
erators were used. The first generated events for a phase- 
space distribution of thep7r°7r° final state (ps). The second 
generator, called 'sequential' (seq), simulated the reaction 
chain 7p Z\(1232)7r° piT^n'^ using a realistic mass 
distribution for the Z\(1232). This was motivated by pre- 
vious results |i22 [ l25 ) . indicating a significant contribution 
from decays of N* resonances to the Z\(1232) intermedi- 
ate state, which gives rise to invariant-mass distributions 
different from phase-space behavior. 

From both simulations the detection efficiency was ex- 
tracted as e{Ery , 771(7^)) and as e(i?^, m(7rp)), i.e. as a 
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Fig. 8. Double tt" simulated detection efficiency for the 7p — )■ 
Ptt^tt" reaction. Main plot (blue) upward triangles: average of 
efficiency extracted from pion-pion and pion-proton invariant 
mass distributions for phase-space event generator, (red) down- 
ward triangles: same for sequential event generator. (Black) 
solid line: average of (ps) and (seq). Insert: deviation of all 
four efficiency curves from the adopted average. 



function of incident photon energy E-y and the invari- 
ant mass m(7r7r) of the pion pairs and the invariant mass 
7Ti(7rp) of the pion-proton pairs. For the latter, the two 
identical pions were randomized. Both types of invariant- 
mass distributions were then corrected with the appro- 
priate detection efficiencies and the detection efficiency 
for the total cross section computed by integration over 
the distributions. The results are summarized in Fig. |51 
Since the efficiencies obtained with the two different event 
generators differ only slightly, their average was used. As 
shown in the insert of Fig. [8] the different results agree 
within ±2% - ±3% at lowest incident photon energies and 
within ±6% at highest incident photon energies. We there- 
fore estimate a systematic uncertainty rising from 3% at 
310 MeV incident photon energy to 6% at 800 MeV. Some 
typical invariant-mass distributions are compared to the 
distributions used in the event generators and to previous 
data in Fig. Uni 

In a more recent measurement with the CB/TAPS 
setup at MAMI-C, which will be published elsewhere |68j, 
data for 27r° production off the proton was also taken up 
to higher incident photon energies (1.4 GeV) and ana- 
lyzed in a different way, using the kinematic-fit technique 
described in [5^155] . Due to this analysis and different 
trigger conditions, the detection efficiency was quite dif- 
ferent from the present experiment (see Fig. [5]) . For the 
range of incident photon energies discussed here, it was al- 
most constant around 60%. The MAMI-C measurements 
did not reach the same statistical precision in the thresh- 
old region as the present results, but they can serve as an 
independent cross-check for systematic uncertainties. 

The efficiency correction is more critical for the 7p — > 
nTT^TT"'' reaction since in this case the charged pions were 



only accepted in the CB, which excludes a small part of the 
reaction phase-space and must be extrapolated by the ef- 
ficiency simulations. For this reaction three different event 
generators were used: phase-space and the reaction chains 

7P ^ Z\°(1232)7r+ ^ n7r°TT+ (8) 
7P ^ Z\+(1232)7r" ^ mr"TT+. (9) 

In a first step, the uncorrected invariant-mass distri- 
butions of the riTr" and the nTr"*" pairs were fitted with a 
superposition of the results from the three different simu- 
lations (the TT^TT^-pairs were much less sensitive to these 
reaction mechanisms). The relative contributions of these 
processes are summarized in Fig. |9l 
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Fig. 9. Relative contribution fd of phase-space (black dots), 
Zi "^-decays (red triangles) and Z\^-decays to the invariant-mass 
distributions of jp — )■ niY^Tv^ . Symbols and curves represent 
the average of the fits for the nvr" and mv^ distributions; the 
error bars indicate their difference. 

It should be emphasized that this analysis was only in- 
tended to construct a realistic detection efficiency. It does 
not include minor contributions like p-meson decays or in- 
terference terms. Nevertheless, the result reflects properly 
the main features of the nir'^'K'^ reaction. Close to 

threshold it is dominated by the 7p Zi*'(1232)7r+ ^> 
nTr'^TT"'" reaction chain. This is due to contributions from 
pion-pole and Z\-Kroll-Rudermann background terms, where 
a Z\''7r"'"-pair is produced at the flrst vertex and the 
subsequently decays into ni:^ . Such diagrams, with the 
two pions interchanged, do not contribute, since the in- 
cident photon couples only to charged pions. The dom- 
inance of 7p — > Z\'^(1232)7r+ rm^Tv^ at low energies 
is clearly visible in the invariant-mass distributions (see 
Fig. EH). At the highest incident photon energies, contri- 
butions from sequential decays of A''* resonances become 
important. Since, due to isospin invariance, the A'^* — > 
Z\°7r+ n7r°7r+ and A^* Z\+7r° 7i7r+7r° decays have 
the same probability, the contributions from and Z\+ 
decays become comparable. 
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Fig. 10. Typical invariant-mass distributions m(7r'', tt") and m(p, tt") for different ranges of incident photon energy for 7p — >■ 
Ptt^tt". (Red) dots: present data, (black) open squares: Wolf et al. [25], solid lines: phase-space, dashed lines: sequential event 
generator. 
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Fig. 11. Typical invariant-mass distributions (red) dots m(7r'', vr"'"), m(n,7r^), m{n,Tv'') for jp — >■ nTr^Tr"*". (Black) squares: 
Langgartner et al. [27]. Solid lines: phase-space, dashed lines: jp Z\''(1232)7r''' — ^ nTr^Tr"*", dotted lines: 7p Z\"*"(1232)7r" — >■ 
mr^-K'^ event generators. 
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In the intermediate energy range, where none of these 
processes is dominant, phase-space behavior parameter- 
izes phenomenologically the contributions from many dif- 
ferent diagrams. The results from the fits of the two dif- 
ferent types of invariant-mass distributions are in quite 
good agreement for most of the energy range; only around 

w 600 MeV do larger deviations occur. 
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Fig. 12. Detection efficiency for the jp — >■ mv^iv'^ reaction. 
(Black) upward triangles: from invariant-mass distributions of 
final state. (Blue) downward triangles: from invariant-mass 
distributions of n'^n final state. (Red) dots: average. Insert: 
deviation from average. (Black) curves: assumed systematic 
uncertainty of efficiency correction. 



For the final efficiency simulations, event generators 
with a corresponding mix of the three contributions were 
used and, as in the 7r°7r° case, the detection efficiency was 
corrected as a function of the photon energy and the 
invariant mass of the particle pairs. As a check for sys- 
tematic effects, this was done independently for the 7r°n 
and 7r"*"n pairs. The total detection efficiency was then 
again computed by integration. The result is summarized 
in Fig. 1121 The efficiencies obtained this way are in good 
agreement above E^ « 700 MeV. They differ by up to 
±10% at intermediate energies and by up to ±20% in the 
threshold region. For the final correction their average was 
used and the solid (black) curves in the insert of Fig. [T^ 
were assumed as systematic uncertainties, i.e. ±5% at 800 
MeV and ±20% at 350 MeV. 

In contrast to double tt*^ production, the detection effi- 
ciency becomes very small in the threshold region because 
the charged pions are absorbed in the target or other mate- 
rial. Therefore it was not possible to analyze this reaction 
close to threshold. 



4 Results and discussion 
4.1 The threshold behavior 

The total cross section for double 7r° production in the 
threshold region is shown in Fig. 1131 The comparison to 
previous results ^511^ demonstrates the enormous progress 
achieved in experimental precision over the last decade. 
While the first measurements of this reaction 21 , 22] could 
not extract any meaningful results for incident photon en- 
ergies below 400 MeV, the data by Wolf et al. |2S] from 
2000 still have statistical uncertainties on the 100% level, 
the results from KotuUa et al. [55] reduced the uncertain- 
ties to the 50% level, and the present results pushed them 
below the 10% level (to the same magnitude as systematic 
uncertainties), allowing for the first time a stringent test 
of model predictions for this reaction. The more recent 
MAMI-C data do not reach the same statistical precision, 
but do not show any systematic deviation from the present 
results. 
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Fig. 13. Total cross section of the jp 
immediate threshold region. (Red) dots: present measurement, 
open (black) triangles [25], open (blue) squares [22, (black 
stars) MAMI-C 68.. Solid line (dotted line): chiral perturba- 
tion theory prediction [2], Eq. ([TJ (dotted line: Eq. |T]) with 
0.9 nb), dashed line: prediction from Valencia model [17], dash- 
dotted: model by Fix and Arenhovel [ID]. Shaded (green) band 
at bottom: systematic uncertainty of present measurement. 



The measured cross sections (Fig. [T5|) are in excellent 
agreement with the prediction from chiral perturbation 
theory [T3|, using the central value for the s-wave cou- 
pling of the Pi 1 (1440) resonance to the double pion chan- 
nel (Eq. ([T]) with 0.6 nb). The calculation in the framework 
of the Valencia model by Roca et al. [T7] somewhat un- 
derestimates the threshold data, whereas the results from 
the model of Fix and Arenhovel [10] (Two-Pion-MAID) 
are much lower. 

Typical invariant-mass distributions for the threshold 
region are summarized in Figs. 1141 and 1151 Very close to 
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Fig. 14. Invariant-mass distributions for 7p ptt^'tv*^ in the 
threshold region. (Red) dots: present measurement, (black) 
open squares KotuUa et al. [28]. Dashed (blue) curves: phase- 
space, solid (red) curves: model of Fix and Arenhovel [40) . 
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Fig. 15. Invariant-mass distribution of the pion pairs from 
the 7p — !■ ptt^tt" reaction. (Red) dots: present measurement. 
Curves show model results from: Bonn-Gatchina model [19| 
(black, dotted). Fix and Arenhovel [40] (red, solid), phase- 
space (renormalized in area, blue, dashed), Valencia model [17] 
(green, long-dashed with FSI, dash-dotted without FSI). Note: 
calculation from [17] is for = 430 MeV. It was re- normalized 
by the ratio of the total cross section at 430 MeV and the av- 
erage of the total cross section between 400 and 430 MeV. 



threshold (see Fig. I14p pion-pion and pion-proton invari- 
ant masses behave like phase-space (blue, dashed curves). 
At slightly higher incident photon energies (see Fig. [15]) 
the pion-pion distribution is still similar to phase-space 
but develops some excess to large invariant masses. This 



is the typical energy range where in-medium modifica- 
tions have been searched for in quasi-free production off 
heavy nuclei .50i.5l] . The results from the reaction models 
are very different for this energy range. The BoGa analy- 
sis [T^ resembles phase-space in shape but overestimates 
the data on an absolute scale. The prediction from the 
Two-Pion-MAID model [40] on the other hand underes- 
timates the data but also has a different shape. The cal- 
culation from Roca et al. [17 is closest to the data and 
also shows an accumulation of strength at high invariant 
masses, although this effect is slightly overestimated. The 
double-hump structure is due to an interference between 
the isospin 7 = and 1 = 2 amplitudes, which is large 
and destructive in the Valencia model [T7] . It is also note- 
worthy that the pion-nucleon final-state interaction in the 
1 = channel has a large effect on the cross section in 
this model. It roughly doubles the result (compare dashed 
and dash-dotted curves in Fig.lTS]) and effectively accounts 
for the loop-corrections in chiral perturbation theory. This 
effect is not included in the Two-Pion-MAID model and 
probably explains at least part of the missing strength in 
this model. 



The threshold behaviors of 7p — ^ 



pn^TT^ and 



are compared in Fig. [THJ Unfortunately, it is not 



possible to measure jp n7r' 7r+ very close to threshold, 
since the low-energy charged pions do not reach the de- 
tector. The only alternative would have been to detect the 
neutral pion and the neutron and identify the reaction via 
overdetermined kinematics, but such events were not in- 
cluded in the trigger conditions (sum of energy deposited 
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Fig. 16. Comparison of the threshold behavior of the jp — > 
Ptt^tt" and jp mv^iv^ reactions. (Red) dots: present 
measurement, (black) stars: MAMI-C for 'yp =>■ ptt^tt". 
(Blue) squares present measurement, open (black) squares 
Langgartner et al. [27] for yp nvr^Tr"'". Dashed curve: Two- 
Pion-MAID model [40] for yp =>■ nTr^vr"'". Shaded (green) band: 
systematic uncertainty of the present yp — >■ nTr^vr"'' measure- 
ment. Insert: comparison on logarithmic scale, solid (black) 
line: ChPT prediction for 7r°7r° [H], dotted (blue) line: ChPT 
prediction for 7rV+ [13]. 
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in Crystal Ball larger than 60 MeV, combined multiplic- 
ity of hits in Crystal Ball and TAPS three or larger). Over 
the range that could be investigated the cross section is 
larger for the mixed-charge channel; however, at the lower 
limit of the accessible range there seems to be some indi- 
cation for the cross-over predicted by chiral perturbation 
theory. A direct comparison to the predictions from chiral 
perturbation theory is not possible because they are lim- 
ited to energies below 335 MeV (see insert of Fig.[Tni). Also 
shown in Fig.[T6]is the prediction from the Two-Pion-Maid 
model |40| . which again underestimates the cross section, 
although not as dramatically as in the tt^tt" case. 



4.2 The resonance region 

4.2.1 Total cross sections and invariant-mass distributions 

The total cross section for double 7r° production is shown 
in Fig. [T7] and the invariant-mass distributions are sum- 
marized in Figs. US] and [501 The agreement between the 
present data and the MAMI-C data [68] is excellent over 
the whole energy range. 
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Fig. 17. Total cross section for the 7P pty^tt" reaction. (Red) 
dots: present measurement, (shaded (green) band: systematic 
uncertainty), (black) stars: MAMI-C [68], (black) squares: Ko- 
tuUa et al. [21], (blue) triangles: Sarantsev et al. [TS], (cyan) 
diamonds: Ahrens et al. ]231- (Black) stars: new MAMI data. 
Solid curve: Valencia model (Nacher et al. [5§), dashed curve: 
Two- Pion- MAID (Fix and Arenhovel [40]), dotted curve: BoGa 
(Sarantsev et al. [19)V 



For the invariant mass distributions only the present 
results are shown, since in most energy bins the MAMI-C 
data fall exactly on top of them. The total cross section 
agrees within the systematic uncertainties with the ear- 
lier measurements with TAPS at MAMI [2g, the GDH 
experiment at MAMI j33j, and for most of the energy 
range (excluding a small region around 550 MeV) also 
with the Crystal Barrel experiment at ELS A For 



energies above 700 MeV, the present results are slightly 
higher than previous measurements, but also still within 
systematic uncertainties. Concerning the systematic ef- 
fects it should be noted that only the present measure- 
ment and the MAMI-C experiment detected the two pi- 
ous over the full reaction phase space with fairly large 
efficiency (30% - 70% for the present measurement (see 
Fig. [5]), w60% in the MAMI-C case) and thus did not need 
any model dependent extrapolations. This is reflected in 
the excellent agreement of the detection efficiencies simu- 
lated with different event generators and applied in differ- 
ent ways (see Fig. [5]). Furthermore, since detection of the 
recoil proton was not required, only the very well under- 
stood Monte Carlo simulation of electromagnetic showers 
induced from photons was needed for the efficiency calcu- 
lation. 

The agreement with the Valencia model li!^ and the 
BoGa analysis ^TU\ is comparable and clearly better than 
with the Two-Pion-MAID model 40^ . The BoGa coupled- 
channel model was fitted (in addition to many other chan- 
nels) also to the previous TAPS- and CBELSA double 
7r° data [12], so that reasonable agreement could be ex- 
pected. Nevertheless, for incident photon energies below 
600 MeV it overestimates the magnitude of the cross sec- 
tion and does not agree well with the shape of the pion- 
pion invariant-mass distributions (Fig. I18p . In the same 
energy range the Two-Pion-MAID model [40] largely un- 
derestimates the magnitude and strongly disagrees with 
the shape of the pion-pion invariant-mass distributions. 
The Valencia model without FSI [TB1I55] is in quite good 
agreement with the data already at incident photon ener- 
gies above 550 MeV (see Fig. [TS]), while below 450 MeV 
the influence of FSI is large (see Fig. [T5)) . All three mod- 
els reproduce pion-pion and pion-proton invariant masses 
quite well at incident photon energies above 700 MeV, i.e. 
in the range where the reaction is supposed to be domi- 
nated by sequential resonance decays via the intermediate 
Z\(1232), the signal of which is clearly visible in the pion- 
proton invariant mass. 

The total cross section for 7p — s> m:^iT^ is shown in 
Fig. [H] and the invariant-mass distributions are summa- 
rized in Figs. [12] and [^l] Agreement between the present 
and previous measurements is within their systematic un- 
certainties. The difference in the peak maximum around 
750 MeV (also within systematic uncertainties) is proba- 
bly due to the different procedures used for the simulation 
of the detection efficiencies. The procedure used in the 
present work is described in Sec. 13.31 Langgartner et al. 
[27] used a simple phase-space model, and Ahrens et al. 
[32] used the Valencia model [T5j . For this channel, agree- 
ment with the results from the Valencia model [35] and 
Two-Pion-MAID [?(T is comparable and for Two-Pion- 
MAID much better than for the double neutral channel. 
An analysis in the framework of the BoGa model is not 
yet available. At low incident photon energies the pion- 
nucleon invariant-mass distributions clearly show an en- 
hancement of the — n7r° decay from background terms 
(which is also reflected in an enhancement of small wk^ 
invariant masses). 
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Fig. 18. Invariant-mass distributions of pion-pion and pion-proton pairs for the 7p — ^ p-K^Tv" reaction for incident photon energies 
from 400-610 MeV. (Black) dots: present measurement, (red) solid curves: Two-Pion-Maid model by Fix and Arenhovel [40) . 
(black) dotted curves: BoGa analysis Sarantsev et al. [TS], (blue) dashed curves: Valencia model Nacher et al. [SS] (only available 
for the highest incident photon energy). 
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Fig. 19. Invariant-mass distributions of pion-pion and pion-neutron pairs for the 7p — >■ nTr^Tr"*" reaction for incident photon 
energies from 400-610 MeV. (Black) dots: present results, (red) solid curves: Two-Pion-Maid model by Fix and Arenhovel |40] . 
(blue) dashed curve: Valencia model Nacher et al. [38] (only available the highest energy bin of pion-pion invariant mass). 
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Fig. 20. Invariant-mass distributions of pion-pion and pion-proton pairs for the 7p — >■ pn n reaction for incident photon energies 
from 610-820 MeV. (Black) dots: present measurement, (red) solid curves: Two-Pion-Maid model by Fix and Arenhovel |40) . 
(black) dotted curves: BoGa analysis Sarantsev et al. [TU], (blue) dashed curves: Valencia model Nacher et al. |38| . 
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Fig. 21. Invariant-mass distributions of pion-pion and pion-neutron pairs for the 7P — ^ n-n^-K^ reaction for incident photon 
energies from 610-820 MeV. (Black) dots: present measurement, (red) solid curves: Two-Pion-Maid model by Fix and Arenhovel 
[40) ■ (blue) dashed curves: Valencia model Nacher et al. [38] (only available for pion-pion invariant mass). 
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At higher incident photon energies the A signal ap- 
pears in both pion-nucleon invariant masses as expected 
for sequential resonance decays. The build-up of strength 
at large values of the pion-pion invariant mass at the high- 
est incident photon energies has been assigned to a con- 
tribution from /9-meson decays |241I27II38IHD] . 
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Fig. 22. Total cross section of the 7p nTr^Tr"*" reac- 
tion. (Blue) squares: present measurement (shaded (green) 
band at bottom: systematic uncertainty), open (black) squares 
Langgartner et al. [271, magenta diamonds: Ahrens et al. 1321 . 
Model results from Valencia (Nacher et al. ^38j) and Two-Pion- 
Maid (Fix and Arenhovel [ID]). 



4.2.2 Beam helicity asymmetries 

Several model predictions [S^HOirTUlITT] indicate that po- 
larization observables are extremely sensitive to the reac- 
tion mechanisms. So far experimental results are scarce. 
Beam asymmetries have been reported for the 7p — > •p'K^'K^ 
P5] and 7n — n7r°7r° [30j reactions from the GRAAL ex- 
periment. The helicity dependence of the cross sections 
has been measured with the Gerasimov-Drell-Hearn project 



at MAMI for the n7r+7r°, ptt^tt^, p7r+7r~, and p7r^7r° final 
states [5^I551[M1[55] . The reaction models compare to all 
these data in a similar way as to the unpolarized data. 
The predictions from the Valencia model [TOIIBB] and the 
Fix and Arenhovel model ID^, without agreeing in detail 
with the data, reproduce the main features of the split 
into (Ti/2 and (T3/2 components. The strong dominance of 
the 173/2 component observed for the -pi^Pifi final state [55] 
disfavors the Laget model [29j where the dominant con- 
tribution is from the Roper resonance. In this situation 
it came as a surprise when the first measurement of the 
beam helicity asymmetry of the 7p — > pvr+Tr^ reaction at 
the CLAS facility [37^ produced results that could not at 
all be reproduced by any reaction model. 

The present experiment has confirmed this result and 
for the first time measured this polarization observable 




Fig. 23. Definition of the angle $ between reaction plane (in- 
coming photon and outgoing nucleon) and production plane 
(pion pair). For identical pions (tt^tt" pairs) the role of pion 
one and two has to be randomized. 



for the 7r°7r° and 7r''7r+ final states. Results have already 
been published in a preceding Letter [36^. Beam helicity 
asymmetries can be measured for reactions with at least 
three particles in the final state with a circularly polar- 
ized photon beam on an unpolarized target. The helicity 
asymmetry /® is defined by: 



1 da^ — da 



da 



da^ 



1 - iV- 



■iv- 



N+ 



(10) 



where da^ is the differential cross section for each of 
the two photon helicity states, and is the degree of 
circular polarization of the photons. For the extraction 
of the asymmetry /® (^, (9^^ , (9^2 , ...) at fixed kinemat- 
ical parameters the cross sections da^ can be replaced 
by the raw count rates N"^ (right hand side of Eq. ITUl) 
since all normalization factors cancel in the ratio. In prin- 
cipal, detection efficiency weighted count rates must be 
used for angle integrated asymmetries. However, due to 
the K, Att coverage of the solid angle, detection efficiencies 
were rather fiat in phase space at fixed incident photon 
energy, so that the effect of the efficiency corrections on 
the asymmetries turned out to be negligible. The photon 
beam was produced from bremsstrahlung of longitudinally 
polarized electrons. In the energy range of interest, polar- 
ization degrees P.y between 60% and 80% were achieved. 
The angle <P between reaction and production plane (see 
Fig. I23p was constructed in the same way as in the work 
of Roca [7T]. For the 7r°7r+ final state, the pions were or- 
dered as shown in the figure, i.e. ^ is the angle between 
the reaction plane and the part of the production plane 
with the charged pion. For the double 7r° final state their 
assignment was randomized, which means that the asym- 
metry must obey /® (<?) — I'^{<P + tt). This was taken into 
account in the modeling, but not enforced in the data 
analysis. 

The advantage of this polarization observable is two- 
fold. It can be measured with good statistical quality (since 
only the electron beam must be polarized) and with small 
systematic uncertainties (since most uncertainties cancel 
in Eq. (jlOp ) and it is very sensitive to different contribu- 
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tions in the reaction models as has been demonstrated in 

The resuhs for the tt'^tt'^ and tt'^tt^ channels are sum- 
marized in Figs. [M] and [25l Parity conservation enforces 
that /®(^) = -/®(27r - (p). This condition was not used 
in the analysis, but as indicated in the figures it is almost 
perfectly respected by the measurements. 
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Fig. 24. Beam helicity asymmetry for 7p — >■ pyr'^Tr'^ for four dif- 
ferent ranges of incident photon energy. Filled symbols: 7®(<?), 
open symbols: — 7®(27r — (Red) solid curves: Two-Pion- 
MAID gO], (blue) dashed curves: Valencia model [IT], (black) 
dotted curves: BoGa fit [IQJ. 



The results for the double 7r° channel are in reasonable 
agreement with the Two-Pion-MAID model [10] and the 
BoGa analysis [12]. The results from the Valencia model 
|17j are completely out of phase with experiment. 

In the case of the 7r''7r+ final state, no model reproduces 
the experimental results. This reaction is more compli- 
cated than 7p pir'^Tr^ , which is dominated in the models 
by sequential resonance decays while jp — nTr^Tr"*" has ad- 
ditional contributions from the p mesons and also stronger 
contributions from non-resonant background terms. The 
predicted infiuence of such terms is shown in Fig. 1251 where 
for the Valencia model the full calculation and the result 
from a truncated model without the contributions from 
p-meson decays are compared. Surprisingly, inclusion of 
the p terms, which substantially improved the agreement 
for the total cross section and the invariant-mass distri- 
butions, has no positive effects for the asymmetry. The 
models do not even come close to the measurement. It is 
interesting to note, that for this reaction also the relative 
contribution of and ai/2 components [32] is quite dif- 
ferent in the Valencia and the Fix and Arenhovel model 
and for both of them agreement with experiment is worse 
than for the other isospin channels. This means that the 
reaction mechanisms are still not understood in detail. 
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on 

Fig. 25. Beam helicity asymmetry for jp — >■ nyr'^Tr^ for six dif- 
ferent ranges of incident photon energy. Notation as in Fig. [21] 
except (black) dash-dotted curves: Valencia model ^1 without 
p-contributions. Results from BoGa are not available. 

5 Summary and Conclusions 

Precise cross-sections and beam helicity asymmetries have 
been extracted for the jp — >■ pn'^Tr^ and 7p — >■ n7r'^7r+ 
reactions from the production thresholds up to the second 
resonance region. 

In the threshold region the results support strongly the 
predictions from chiral perturbation theory [14| . The to- 
tal cross section for double 7r° production agrees with the 
ChPT prediction within statistical uncertainties. It sup- 
ports also the value used in the ChPT calculation for the 
s-wave coupling of the Pii(1440) resonance to the dou- 
ble TT channel. The cross section calculated in the frame- 
work of the Valencia model [TTj, taking into account wN 
re-scattering effects that resemble the loop corrections of 
ChPT, is also close to the measurement but underesti- 
mates it slightly outside the experimental uncertainties. 
On the other hand, the effective Lagrangian model (Two- 
Pion-MAID) from Fix and Arenhovel [ID] underestimates 
the threshold cross section typically by a factor of five. 
This underlines the large importance of ttA'^ re-scattering 
for the double tt*^ channel. 

The results for 7r"7r+ could not be directly compared 
to the chiral perturbation prediction. There is a gap of 
«15 MeV between the upper energy limit of the theory 
prediction and the lower energy limit of the experimental 
sensitivity. Nevertheless, below incident photon energies 
of 400 MeV the 7r°7r+ excitation function bends down- 
wards, seems to cross the tt^tt" data around 350 MeV, and 
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could approach the ChPT threshold prediction. But it is 
of course desirable to close the gap, which in terms of cross 
section spans almost two orders of magnitude, either by 
more refined ChPT calculations reaching higher incident 
photon energies or by more sensitive experiments. Alto- 
gether, the experimental results clearly support the ChPT 
calculations four threshold production of pion pairs. 

The situation is much less clear at higher incident pho- 
ton energies where the contributions from nucleon reso- 
nance decays become important. Fairly large discrepancies 
between experiment and model results and between differ- 
ent models occur in the intermediate energy range from 
400-650 MeV, where contributions from the Pii(1440) res- 
onance have been discussed [19]. In the tt^tt" channel the 
Two-Pion-MAID model 'W' strongly underestimates the 
magnitude of the cross section while the BoGa analysis 
[T9] overestimates it. Both models disagree also in dif- 
ferent ways with the shape of the invariant-mass distri- 
butions. Agreement for Two-Pion-MAID with the 7r°7r+ 
results is better as far as the magnitude is concerned but 
also here the shape of the invariant-mass distributions is 
not well reproduced. The results for total cross sections 
and invariant-mass distributions from the Valencia model 
are at least in reasonable agreement with experiment over 
the whole energy range. In the version including nN re- 
scattering, the tt'^tt" invariant-mass distributions are quite 
well reproduced even at very low incident photon energies. 

Finally, at the highest incident photon energies, in the 
second resonance region, all model analyses reproduce the 
absolute magnitude and the main features of the shape 
of the invariant-mass distributions for both isospin chan- 
nels. However, the models still do not agree on the rela- 
tive importance of the different contributions, e.g. for the 
-033(1700) resonance to double 7r°. The contribution of 
this state to tt^tt" is very strong in the BoGa analysis, 
much smaller in Two-Pion-MAID, and almost negligible 
in the Valencia model. 

As a new tool, beam helicity asymmetries, which had 
been predicted to be very sensitive to the details of the 
models, [TTj have been measured with high precision. The 
result is surprising. The Valencia model, which had the 
best overall agreement with the measured cross section, 
could not reproduce this observable for any of the isospin 
channels. Two-Pion-MAID and the BoGa analysis did much 
better for the double 7r° channel, but so far no model could 
reproduce the results for the Tr^Tr"*" channel (an analysis in 
the framework of BoGa is not available for this channel). 
Therefore we must conclude, that so far none of the avail- 
able reaction models correctly reflects the details of double 
pion photoproduction in the resonance region. Further ef- 
forts on the theory side are highly desirable. 
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